Molecular dynamics study of hcp/fcc nucleation and growth in bcc iron driven by 

uniaxial strain 
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Molecular dynamics simulations are performed to investigate the structural phase transition in 
body-centered cubic (bcc) single crystal iron under high strain rate loading. We study the nucleation 
and growth of the hexagonal-close-packed (hep) and face-centered-cubic (fee) phases, and their 
crystal orientation dependence. Results reveal that the transition pressures are less dependent on 
the crystal orientations (^^14 GPa for loading along [001], [Oil], and [111] directions). However, 
the pressure interval of mixed phase for [Oil] loading is much shorter than loading along other 
orientations. And the temperature increased amplitude for [001] loading is evidently lower than 
other orientations. The hcp/fcc nucleation process is presented by the topological medium-range- 
order analysis. For loading along [001] direction, we find that the hep structure occurs firstly 
and grows into laminar morphology in the (Oll)bcc planes with a little fee atoms as intermediate 
structure. For loading along [Oil] and [111] directions, both the hep and fee structures nucleation 
and growth along the {llOjbcc planes are observed, whose morphology is also discussed. 

PACS numbers: 64.70.kd, 71.15.Pd, 61.50.Ks, 62.50.-p 
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INTRODUCTION 



Besides melting, dislocation generation, and twinning, 
structural phase transition is also one important physical 
process in understanding the plasticity of solid (metallic, 
particularly) materials. Among various metallic mate- 
rials, iron, due to its critical importance in condensed- 
matter physics, materials science, geophysics, and human 
development, have been extensively studied in aiming at 
understanding its structures and phase transition under 
extreme conditions. Thereinto, the research on classi- 
cal phase transition of iron from ferromagnetic body- 
centered cubic (bcc) ground state a phase to nonmag- 
netic hexagonal-close-packed (hep) e phase has a history 
of more than half a century. Early representative work 
was carried out in 1956 by Bancroft et al. [1], where 
they discovered the transition under shock wave com- 
pression at pressure around 13 GPa. Thereafter, plenty 
of theoretical and experimental works have been done to 
study the structure properties, inner physics, and me- 
chanical features involving this a-^e transition process 
0, 11 II i, H, 0, Hi- Among these studies, in 1962, 
Jamieson and Lawson obtained crystallographic evidence 
for the a— >e transition at 13 GPa in static x-ray diffrac- 
tion experiments [2,], which confirmed previous shock- 
wave experiment of Bancroft et al. [ij; in 1991, Tay- 
lor and Pasternak reported the onset of the transition 
in range of 9-15 GPa with a large hysteresis width of 
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about 6 GPa [4|; in 1997, Boettger and Wallace presented 
theoretical analysis of metastability and dynamics of the 
transition in 1999, Wang and Ingalls suggested three 
possible models for the a— >e transition mechanism 0]; 
and in 2004, Caspersen et al. theoretically investigated 
the influence of shear in the transition based on a multi- 
scale model [1]. 

However, the structural transition is a mutation pro- 
cess, which leads to direct observation of the phase transi- 
tion very difficult in experiment. And the abrupt changes 
in physical and mechanical characteristics also bring 
about difficulty in theoretically analyzing the dynamic 
response and micro-behavior of materials. Therefore, it 
is necessary to carry out systematic simulations at atomic 
level by using molecular dynamics (MD). MD is a method 
for computationally evaluating the thermodynamic and 
transport properties of materials by solving the classical 
equations of motion at the atomic level. Many preceding 
works have testified that MD is the most conceptually 
straightforward method and an altogether very utilitar- 
ian method for obtaining atomic details. In 2002, Kadau 
et al. [9] realized the micro-simulation of shock-induced 
phase transition of solid iron along [001] direction by em- 
ploying MD method with an embedded atom method 
(EAM) potential [l3, El , showing the micro- features of 
a— >e transition. Subsequent ultrafast x -ray diffraction 
experiments processed by Kalantar et al. 11211 v erified the 
validity of this work. In 2005, Kadau et al. (ISl] simulated 
shock-induced phase transformations in bcc iron induced 
by shock loading along [001], [Oil], and [111] directions 
to research the orientation dependence of the developing 
microstructure on the crystallographic shock direction. 
They found that the fee phase appeared especially when 
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loading_along [Oil] and [111] directions. In 2007, Kadau 
et al. [ij] utilized the MD method again to simulate the 
polycrystalline iron in shock loading condition, simula- 
tion results correspond properly with experiment ^15<] . 

These abovementioned MD studies have greatly en- 
riched our knowledge of the structural transition of iron 
and, needless to say, subsequently inspired more simu- 
lation works [l^, [13] ■ In our last work we have 
simulated the a—>e structural transition of iron under 
isothermal compression along the [001] direction. It was 
showed that the laminar structure forms along {llOjbcc 
planes. On the other hand, when the loading direction 
is different, both the hep and fee nucleation and growth 
in bcc single crystal iron will occur. To date, no detailed 
investigation has been proceeded to study the orientation 
dependence of hep and fee nucleation and growth. This 
issue is focus of our present paper. Specially, by adia- 
batic compression along [001], [Oil], and [111] directions, 
in this paper, we investigate the a— *-e/7 structural phase 
transition in bcc single crystal iron. The pressure and 
temperature effect, nucleation mechanism and morphol- 
ogy characters are all discussed. 

The rest of this paper is arranged as follows. In Sec. II 
the simulation method is briefly described. In Sec. Ill we 
present and discuss our simulated results. A summary of 
our main results is given in Sec. IV. 



II. METHODOLOGY OF THE SIMULATIONS 

For this work, we use classical MD and the EAM 
potential of iron developed by Voter-Chen (VC) [l^ 
to conduct the simulations. Our calculations are per- 
formed in three cases by using three samples A, B, 
and C to simulate the loadings along the [001], [Oil], 
and [111], respectively. Specially, sample A consists of 
SOj. [100] x50j, [010] ^100^ [001] atomic cells, sample B con- 



cells, and sample C 



[111] 



cells. Periodical 



sists of 36j;,[oii]XlOOj,jioo] x722_[oii] 
consists of 362; [110] x40j^ [JI2] xl80^ 

boundary conditions are adopted in all directions. The 
original temperature in all samples is set at 60 K by speed 
calibration method [l^ . The initial state is manipulated 
at zero pressure through modulating the lattice constant. 
Here the lattice constant of all samples is ao=0. 28725 nm. 
Loading with high strain rate of ~10® along the [001], 
[Oil], and [111] directions is implemented by shortening 
the lattice constants of these three orientations step by 
step during the nonequilibrium MD process. The Ver- 
let algorithm 2Ql is used to integrate the equation of 
motion. The simulation time is up to 200 ps, with the 
time step set at 2 fs. Special techniques are required to 
analyze the huge amount of data produced in the MD 
simulations. Here, the local structure around each atom 
are resolved by using the topological medium-range-order 
analysis [2l[ technique. The phase mass fraction analy- 
sis and radial distribution function are also presented to 
obtain the phase transition information. 
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FIG. 1; (Color online) Variation of pressure with the strain e 
for loading along [001], [Oil], and [111] directions. The mixed 
phase intervals start at Tajooi] (0.072, 14.21), Ta_[oii] (0.071, 
14.66), and T.jni] (0.066, 14.16) and finish at Tf^poi] (0.104, 
19.25), Tf,[oii] (0.075, 13.17), and Tf,[iii] (0.088, ' 16.74) for 
[001], [Oil], and [111] loadings, respectively. 




FIG. 2; (Color online) Variation of temperature with the 
strain e for loading along [001], [Oil], and [111] directions. 
The letters Ts,i and Tt,i (? stands for [001], [Oil], and [111]) 
have same strains with that in Fig. [1] 



III. RESULTS AND DISCUSSION 

A. Pressure and temperature effects 

To begin with, we perform the pressure analysis dur- 
ing the loading process. The calculated pressure as a 
function of compression strain e is shown in Fig. [U 
where e=l-V/Vo and Vq is the initial volume of the 
samples. With increasing the compression strain, the 
following prominent features can be seen from Fig. [1] 
(i) The transition critical points in the strain-pressure 
space are T.jooi] (0.072, 14.21), T^^^oii] (0.071, 14.66), and 
Ts,[iii] (0.066, 14.16) for loading along [001], [Oil], and 
[111] directions, respectively, which illustrates that the 
critical pressures of transition are less dependent on the 
crystal orientations. This is similar to the recent shock- 
loading simulations [13]. However, the strain of Tg [mj is 
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FIG. 3: (Color online) Mass fractions A for the fee and hep 
phases change with strain e for loading along [001] (a), [Oil] 
(b), and [111] (c) directions. The letters Ta,i and Tf_i {i stands 
for [001], [Oil], and [111]) have same strains with that in Fig. 
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evidently smaller than that of other two critical points, 
suggesting that sample C begins transition first. After 
the critical points the pressure is lowered due to the nu- 
cleation of the new phases; (ii) Close to the transition 
critical point Tsjm], the pressure-strain curve for [001] 
loading is below other two curves. This indicates that 
samples B and C experience a more evident hardening 
process compared with sample A before phase transi- 
tion; (iii) The mixed phase intervals are Tg [qqi] to Tf [ooi], 
Tsjoii] to Tfjoii], and Ts,[iii] to Tf^m] for loading along 
[001], [Oil], and [111] directions, respectively. Obviously, 
the mixed phase interval in [Oil] loading is much shorter 
compared with that in [001] and [111] loadings. This is 
due to the fact that the compression direction in sam- 
ple B is parallel to the nucleation planes (see below for 
further discussion). 

Consistent with previous reports [l^, [12] , in our sim- 
ulations we have also observed two main steps during 
structural transition of iron: firstly, atoms in nucleation 
planes {Olljbcc are compressed to form a hexagonal 
structure; secondly, neighbor {110}bcc layers experience 
a relative slip to transform into hep structure. 



Under adiabatic compression, the temperature is set 
free, i.e., not controlled artificially. Thus the tempera- 
ture evolution with the strain can be calculated based on 
our data. Results are presented in Fig. [51 One can see 
that the temperature values at the three starting transi- 
tion critical points (namely, Tg [qqij, Tg ^on], and Tg [mj) 
are all near 64 K, which further implies that the onset of 
phase transition is insensitive to the loading orientations. 
Right after the onset of the phase transition, however, 
loadings along the three different directions display re- 
markably different temperature effects by the fact shown 
in Fig. [2] that the temperature is higher one by one with 
a sequence of [001], [Oil], and [111] in loading orienta- 
tion. While the temperatures of samples B and C change 
up to 185 K and 223 K, the temperature of sample A 
only increases to 102 K at the end of loading. Also, the 
temperature increase style is different. The temperature 
of sample A increases in an up-stair-like manner, which 
can be seen more carefully near the labels Tgjooi] (0.072, 
63.63), Tfjooi] (0.104, 81.17), and followed by a sawtooth- 
like slope after exceeding the strain of 0.136. Here, the 
temperature climb up the first stair may contribute to the 
first step of transition mechanism. This up-stair-like pro- 
cess only costs about 2 ps. Then, the temperature keeps 
almost no change until the second stair. The second and 
third stairs in the temperature of sample A have strong 
relation with the second step in the structural transition 
of iron. Closer view can find that the temperature fluc- 
tuates dramatically after the third stair and only mildly 
fluctuates after the first two stairs. These features of 
temperature in sample A has not been observed in sam- 
ples B and C. Instead, the temperatures in samples B 
and C first alter suddenly up to 110 K and 142 K from 
the initial value of ^64 K, respectively in 5.3 and 9.3 ps. 
Then they increase steadily with strain. On the whole, 
the remarkable distinction in temperature before and af- 
ter phase transition illustrates the substantial structural 
and energetic complexity during transition process. 

Although the main transition mechanism remains the 
same for the three kinds of loadings, the details are signif- 
icantly different. This can be seen from the phase mass 
fraction analysis. Figures 3(a)-(c) show the mass fraction 
A for the hep and fee phases during their evolution with 
the loading strain along the three directions, respectively. 
In the overall view, the variations of the mass fraction of 
close packed (fcc/hcp) phases as a function of strain in all 
three samples are similar to that of temperature. There- 
fore, the temperature and the nucleation of close packed 
phases have tight relation. In addition, while most atoms 
in sample A nucleate into the hep structure, in samples B 
and C there will emerge near half proportion of fee atoms. 
This is similar to the observation in shock-wave loading 
simulation ^14], where it was found that the hcp/fcc ra- 
tio within a grain decreases the more the shock direction 
deviates from the [001]bcc direction of the initial poly- 
crystal. In our study, with the pressure increasing, more 
fee atoms will form when loading along [Oil] and [111] 
directions. The hep mass fraction decreases with the in- 



FIG. 4: (Color online) Morphology evolution of the fee phase 
(the first row) and the hep phase (the second row) in sample 
A with strains labeled under each snapshot. The coordinates 
are shown in panel (a) and the compression directions is along 
Z axis. Color coding was obtained by using the topological 
medium-range-order analysis; white: fee, red: hep, yellow: 
grain boundaries. 




FIG. 5: (Color online) Morphology evolution of the fee phase 
(the first row) and the hep phase (the second row) in sample 
B with strains labeled under each snapshot. The coordinates 
are shown in panel (a) and the compression directions is along 
X aM.s. Color criterion same with Fig. 

crease of fee phase mass fraetion, which indicates that 
the nucleation of fee is due to the slip of neighboring 
hexagonal layers. 

B. Nucleation and growth of hep and fee phases 

The nucleation and growth of the hep phase under [001] 
compression has been investigated in Ref. Here, 
we focus on the emergence and growth of both hep and 
fee phases along different crystal directions by continu- 
ous uniaxial compression. Over the critical pressure, the 
close packed phases begin to nucleate in all three kinds 
of loading. The microscopic view of the initial nucleation 
and growth of the fee (upper panels) and hep (lower pan- 
nels) phases for samples A, B, and C are shown in Figs. 
4-6, respectively. The evolved loading strain correspond- 
ing to these snapshots are also labeled. For an over- 
all view, the morphologies of samples A, B, and C are 
also shown in Fig. [7] at some evolved strains. Based on 
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FIG. 6: (Color online) Morphology evolution of the fee phase 
(the first row) and the hep phase (the second row) in sample 
C with strains labeled under each snapshot. The coordinates 
are shown in panel (a) and the compression directions is along 
Z axis. Color criterion same with Fig. |4] 
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FIG. 7: (Color online) Compressed samples under different 
strains with the first row (three snapshots) for [001] loading, 
the second row for [Oil] loading, and the third row for [111] 
loading. The compression directions of all samples are from 
top to bottom. The coordinates are shown in (a), (d) and 
(g) for loading along [001], [Oil], and [111] directions, respec- 
tively. The green atoms stand for the bcc structure and other 
color criterion same with Fig. [l] 



these figures, we can clearly see that for the [001] and 
[111] loadings, until strain being loaded up to 0.09, the 
bcc structure atoms still maintain bulk morphology, as 
shown in Fig. [TUb) and 7(g). For the loading along [Oil] 
direction, on the contrary, only a few bcc atoms disper- 
sively distribute among bulk hep and fee atoms when the 
loaded strain is over 0.075. This further implies that the 
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mixed phase interval in sample B is shorter than in other 
two samples. 

For the [001] loading, as shown in Fig. [H some fee- 
phase atoms only appear in mixed phase interval and 
nucleate on (llO)bcc plane. With increasing pressure, 
the fee atoms will transform to hep or grain boundary 
atoms. From Fig. IHd), we can find that a little amount 
of fee atoms appear only in the twin boundaries. For hep 
atoms, initially, they homogeneously nucleate [Fig. Hl^e)], 
then evolve into irregular grains [Fig. (H^f)]. These grains 
subsequently develops into the typical laminar structure 
along the (Oll)bcc planes, as shown in Fig. HJ^g). Ul- 
timately, the laminar structure grows into larger and 
reugular grains divided by the grain boundaries [Fig. 
IDJh)]. The grain boundary planes are (lOO)bcc planes, 
paralleling with the compression direction. 

While loading along [Oil] (see Fig. [5]), the fee phase 
nucleates on (llO)bcc, (lOl)bcc, and (liO)bcc planes and 
the hep on (lOl)bcc and (liO)bcc planes. These nucle- 
ation planes all belong to the {llOjbcc family. More spe- 
cially, the fee atoms first nucleate as a flaky structure 
[Fig. [3a)], then, together with hep atoms, evolve into 
the thin samdwich structure perpendicular to the com- 
pression direction [see Fig. [S{b),[5lc),[5](f),[5jg), and Fig. 
[TJd)]. At the end of loading, the fee and hep atoms form 
stacking fault structure and relatively thick samdwich 
structure [Fig. O^d) andl^Ih)]. Therein, the samdwich 
structure is divided by the grain boundaries, which can 
be more clearly seen from the yellow atoms in Fig. [T] 
Note that the grain boundary planes are (Oll)bcc planes 
and perpendicular to the compression direction. 

For loading along [111] (see Fig. the nueleation 
and growth of the close packed phases are different. Just 
over the critical pressure, both fee and hep atoms nucle- 
ate to form lots of small but long flaky structure [Fig. 
mja), EKb), [nUe), andHKf)]. From Fig. EKc), a relaxation 
process right after the critical point Tg [mj has been ob- 
served. Although the fractions of fee and hep phases 
are reduced, the arrangement of atoms becomes more 
inerratic, as shown in Fig. ^c) and[6Kg). At the end 
of loading, the stacking fault arrangement comes out in 
the interior of thick samdwich structure. Compared to 
the [Oil] loading, both fee and hep nueleation crystalline 
planes in [111] loading are (lOl)bcc and (Oll)bcc- Besides, 
the grain boundary planes (IlO)bcc m [111] loading are 
parallel, instead of perpendicular (in the case of [Oil] 
loading), to the compression orientation. 

On the whole, for [001] loading, the a^e transition 
happens in the whole body of the sample and fee struc- 
ture only appears as the intermediate structure. How- 
ever, for [Oil] and [111] loadings, the fee and hep struc- 
tures nucleate simultaneously. Note that for [Oil] loading 
the nueleation is localized. The hep structure reaches its 
maximum at the evolved strain of '^0.09 and '^0.1 for 
[Oil] and [111] loadings, respectively. Further increase 
of strain will result in the decease (increase) of hep (fee) 
weight. A visible difference between [Oil] and [111] load- 
ings is that for the latter, there occurs an evident relax- 
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FIG. 8: Radial distribution functions (RDFs) under different 
strains e for crystals compressed along [001] (a), [Oil] (b), and 
[111] (c) directions. As a reference the ideal positions of the 
fee and hep structures, as well as the distribution function for 
the unstrained bcc (e=0), are shown. 



ation process after the critical point. 



C. Radial distribution function analysis 

The slip investigation of atoms must presume upon ra- 
dial distribution function (RDF) analysis and phase mass 
fraction analysis. Figure [S] displays the RDFs under dif- 
ferent compressions, where the strains of all three sam- 
ples are from to 0.15. The ideal positions of fee and 
hep structures have also been plotted. Along with the 
increase of strain the first peak in all three samples devi- 
ates towards left gradually, while the second peak devi- 
ates towards the first peak and disappears finally. After 
the strain becomes bigger than 0.0975, along with the 
disappearance of the second peak, a new peak appears, 
with the position (shown by the arrows in Fig. [8]) ex- 
actly the same as that of the second peak of the fee/hep 
phases. Therefore, these changes in RDFs definitely sig- 
nify atomic transformation to the close packed phases. 
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Moreover, the evident broadening of the peaks in [Oil] 
and [111] loadings after transition is caused by the tem- 
perature effect. The third and fourth peaks at £=0.15 
in Fig. [Hl^a) obviously trend to split; this is prominently 
different from that in Figs. [DJb) and[8ljc). This further 
illustrates that sample A is mainly featured by the hep 
phase at the end of loading. All these analyses indicate 
that the phase distribution and the mechanism of local 
slip in three samples are different in the details of tran- 
sition process. 

IV. CONCLUSIONS 

In summary, we have simulated the process of a— s-e/7 
structural phase transition in bcc single crystal iron un- 
der uniaxial compression along three typical bcc orienta- 
tions by MD method. It is found that samples B and C 
experience a more evident hardening process compared 
with sample A before phase transition. The tempera- 
tures for [Oil] and [111] loadings increase more highly 
than that for [001] loading. The mixed phase interval of 
sample B is much shorter than other two samples, which 
is caused by the fact that the compression orientation be- 
longs to the nucleation planes {110}bcc- The transition 



pressure is the same for all three samples with a value of 
~14 GPa, over which the nucleation of hep and fee atoms 
appears and grows into different morphologies. For [001] 
loading, the a— >e transition happens in the whole body 
of the sample and hep atoms nucleate on (Oll)bcc planes 
to form laminar structure; the fee phase only appears in- 
termediately. For [Oil] and [111] loadings, however, (i) 
nucleation of fee and hep structures occurs simultane- 
ously; (ii) the hep mass fraction reaches its maximum at 
the evolved strain of ^0.09 and ~0.1 and further increase 
of strain will result in the decease (increase) of hep (fee) 
weight; (iii) the fee and hep atoms form stacking fault 
structure and relatively thick samdwich structure at the 
end of loading. Note that for [Oil] loading the nucleation 
is localized. The phase distribution and the mechanism 
of local slip in three samples are different in detail. 
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